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Abstract. In de Sitter spacetime there exists an absolute minimum for the mass of a 
spin-2 field set by the Higuchi bound m 2 > 2H 2 . We generalize this bound to arbitrary 
spatially flat FRW geometries in the context of the recently proposed ghost-free mod- 
els of Massive Gravity with an FRW reference metric, by performing a Hamiltonian 
analysis for cosmological perturbations. We find that the bound generically indicates 
that spatially flat FRW solutions in FRW massive gravity, which exhibit a Vainshtein 
mechanism in the background as required by consistency with observations, imply that 
the helicity zero mode is a ghost. In contradistinction to previous works, the tension 
between the Higuchi bound and the Vainshtein mechanism is equally strong regardless 
of the equation of state for matter. 
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1 Introduction 

Theories that attempt to modify General Relativity (GR) have a long and rich his- 
tory. Their study in most recent years has been further motivated by the observational 
finding on the Supernova data [1, 2] which points to acceleration in the current expan- 
sion of the Universe. If GR is correct, then there must be a dark energy density of 
p ~ 10 _29 g/cm 3 . If this value is due to the cosmological constant A, it will enforce an 
extremely small number for A, very far from the value that arises from quantum field 
theory considerations [3]. On the other hand, if one is willing to modify GR, it has 
been shown in different modified gravity theories [4]- [7] that there exist self-accelerating 
solutions [8, 9] which might explain the current acceleration of the Universe without 
resorting to dark energy. A ubiquitous issue in theories of gravity, such as the smallness 
of the cosmological constant, can have different manifestations in different theories. A 
much desired feature of such theories would be a mechanism by which the smallness of 
the parameter A (or whatever parameter takes its place) is technically natural. Quite 
generally, if setting to zero a small parameter in a theory results in an additional sym- 
metry, it is reasonable to expect that quantum corrections to that parameter will be of 
the same order of the parameter itself as they are protected by the initial symmetry; 
if this is the case, the parameter is said technically natural. 

Take for example massive gravity (MG): at large distances MG weakens with 
respect to GR, the potential reads ~ e~ mr 'jr. This results in the possibility of an 
accelerated expansion without dark energy. Observations force one to assume a very 
small m and again, the fine tuning issue reappears in the ratio m/M-p\. Interestingly 
though, one notices here how the m = theory, GR, has an additional gauge symmetry 
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that is broken by the massive term and so there is indeed room for a technically natural 
small value of the parameter m. 

Departing from GR comes with a heavy baggage: since the pioneering proposal 
of Fierz and Pauli [10], theories of massive gravity have been plagued with continuity 
issues [11-13], instabilities (ghosts) [14], and with a hierarchy of scales which has made 
it hard to make sense of the theory [15]. A lot has been done in recent years to make the 
prospect of a massive theory of gravity more intriguing and, possibly, more predictive. 
This effort [16] culminated [18, 20] in the formulation of a theory of massive gravity 
(which we shall refer to as dRGT from now on) which is free of ghosts at the fully 
non linear level [21-28]. This theory is endowed with a benevolent hierarchy of scales 
which neatly splits the linear regime, the non-linear one and the regime where quantum 
effects must be taken into account. 

Having such a young theory at one's disposal, there is no shortage of aspects 
in need of careful investigation [29-36]. The more cosmologically inclined might for 
example opt for a study of realistic cosmological solutions for dRGT, this has recently 
been done in [37]. The analysis we present here shares some features with the work [37] 
in that it points towards the same direction. In this manuscript we study the classical 
stability of the scalar sector of dRGT theory up to quadratic order in perturbations 
when the reference metric is taken to be FRW. We probe several cases according to the 
background value of the dynamical metric g^ u and the reference metric f^ u . Matter 
content besides the cosmological constant is also considered. Such a study has been 
performed in the past for theories which were known to have ghosts [38, 39]. In this 
context the so called Higuchi bound [40] on the mass of the graviton was introduced. 
It is a strong lower bound on m, m 2 > 2H 2 , and arises from the requirement that the 
kinetic term for the helicity zero mode, i.e. the scalar cosmological perturbations, is 
positive definite. Things do not necessarily improve when one includes matter because 
a similar bound must hold over the different cosmological epochs [39], as we will see. 
We show that, even if in dRGT there is more room to accommodate for such a bound 
(e.g. dRGT has 2 free parameters) and even if we employ at full the freedom on the 
reference metric, once observations are taken into account the bound remains quite 
stringent. We hint to a possible resolution in the Conclusions. Our final bound is give 
by 

H ( H H 2 \ 

m 2 (H) = m 2 — (3 + 3a 3 + a 4 ) - 2(1 + 2a 3 + a 4 )— + (a 3 + 04)772 > 2r2 - 
-no V h o tl J 



where H is the Hubble rate of the dynamical metric and Hq that of the reference (non- 
dynamical) metric, and a 3 and a 4 are the two free parameters in the dRGT model 
[20]. This should be compared with the associated Friedmann equation 

1 ui 2 H 

H 2 = ^T72 P ~ ( 6 + 4 «3 + "4)— + (3 + 3a 3 + a 4 )m 2 — - 

x 2H 2 , ,m 2 H 3 
(1 + 2a 3 + a A )m — ^ + (a 3 + a 4 ) — —3. 

li 6 tt 
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In this paper the background metrics studied are both homogeneous and isotropic, 
a natural next step would then be to consider introducing inhomogeneities at the 
background level. There is more to support further work in this direction: in [37], it 
was shown that there exist no truly homogeneous and isotropic cosmological solutions 
in dRGT theory with a Minkowski reference metric, for the very same reasons that 
guarantee no Boulware-Deser ghost is present in such a theory. On the other hand, 
there are approximate solutions that well describe observations and yet evade this 
no-go theorem. We will expand upon this in a forthcoming work [45]. 

This paper is organized as follows: in the first Section we briefly introduce dRGT 
theory. In Section 2 we report the details of the theory at second order in perturbation 
for the scalar sector. In Section 3 we briefly introduce the Higuchi bound and summa- 
rize previous work on the subject. In Section 4 we present the analysis for the theory 
when both metrics are de Sitter. In Section 5 we add matter content and consider 
FRW solutions. In the Conclusions section we elaborate on our findings and future 
work. 

Although in this work we concentrate on quite specific aspects of the theory, it is 
important to appreciate that this model of massive gravity seems to enjoy properties 
and structures that are ubiquitous in current analysis of, for example, alternative 
models to inflation [47-51]. This theory as a whole also appears to be part of a larger 
family of massive theories of gravity [52] some of which first emerged in the study of 
AdS3/CFT 2 correspondence. 

2 dRGT Massive Gravity 

The theory of massive gravity defined on an arbitrary reference metric f^ u [53] is just 
a straightforward generalization of the theory proposed in [20]. The Lagrangian takes 
the form of Einstein gravity with matter plus a potential that is a scalar function of 
the two metrics 

C = ^f^g ( ( % + 2m 2 U(g, /)) + C M . (2.1) 

The most general potential U that has no ghosts is build out of characteristic polyno- 
mials of the eigenvalues of the tensor 

W(9J) = %-Vfi i ^, (2-2) 

so that 

U(g,H)=U 2 + azUz + c u U A , (2.3) 
where the a n are free parameters, and 

U2 = \ ([IC] 2 - [/C 2 ]) , (2.4) 
W 3 = ^([/C] 3 -3[/C][/C 2 ]+2[/C 3 ]), (2.5) 
U A = i ([/C] 4 - 6[/C 2 ][/C] 2 + 8[/C 3 ][/C] + 3[/C 2 ] 2 - 6[/C 4 ]) , (2.6) 
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where [. . .] represents the trace of a tensor with respect to the metric g^. The absence 
of ghost for this theory for a Minkowski background metric was shown in the decoupling 
limit in [17, 18, 20], fully non-linearly beyond the decoupling limit in [21, 22], as well 
as in the Stiickelberg and helicity languages in [23-28]. Varying with respect to the 
metric g^ u we find the equations of motion 

Gfiu + m 2 X flv = Mpf , (2.7) 

where 

X, v = K» v - K 9lw - (1 + a 3 ) - KK^ + l -g, v ([/C] 2 - [/C 2 ])) (2.8) 

+(a 3 + a*) Ulu - KK?^ + \lC, v ([/C] 2 - [/C 2 ]) - l 9lw ([/C] 3 - 3[/C][/C 2 ] + 2[/C 3 ])) . 

Using the Bianchi identities, we obtain the following constraint on the metric 

m 2 V^X^ = 0. (2.9) 

It is clear from the above formula that the subspaces of parameter choices a- 3 + = 
and o?4 = = 1 are special. We shall see below that the same structure arises in the 
Higuchi bound. As is well understood, we can introduce Stiickelberg fields to recover 
diffeomorphism invariance, however in what follows we shall work entirely in unitary 
gauge to make clearer the comparison with the work of [38] and [39]. 



3 Action at 2nd order in perturbations 

Our starting point for the two metrics reads: 

B _(-N 2 + N?/a 2 Nj \ _(-M\t) 0, \ 

where, being the fixed, reference metric, is not perturbed 1 . It is convenient to em- 
ploy throughout this work the ADM decomposition for the metric and its fluctuations 
[54]. For future convenience we define the ratio of scale factor 

b , x 

- = l + z. 3.2 
a 

If both the background and reference metrics are de Sitter then we would have g^ v = 
(1 + z)f liv with z time independent. However our formula will be valid in the general 
time-dependent case. 

The same constraint, coming from the consistency of the Friedmann a and acceler- 
ation a equations that forbids the existence of FRW solutions for Minkowski reference 
metric, in the generic FRW case imposes the condition 

(1 + 2(1 + a 3 )z + (a 3 + a 4 )z 2 ) ( - - J- ) = 0. (3.3) 

\a H j 



1 Giving / full dynamics would take us into the realms of multi- metric ghost-free theories [58, 59]. 
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where Hq is the Hubble rate for the reference metric Hq = b/(Mb). Although it appears 
that we could solve this relation with (1 + 2(1 + a 3 )z + (a 3 + a^z 2 ) = 0, it turns out 
that in this branch of solutions there are fewer than the full 5 propagating solutions in 
the gravity sector, and these correspond to singular points in phase space which would 
be infinitely strongly coupled due to the vanishing of the kinetic term for the helicity 
zero mode at quadratic order (we shall discuss this further in [44]). 

The correct solution of this relation that gives rise to the full 5 propagating degrees 
of freedom in the gravity sector is 

b H , . 

and it is this solution only that we shall consider in the following. This relationship is 
crucially important in what follows. It is this relation that ties together the dynamics 
of the two metrics. This same relationship arises in the bigravity case when both 
metrics are taken to be dynamical [59] (so does the same issue with the pathological 
branch) . 

In what follows we shall also set Mp\ = 1 to simplify the calculation, reintroducing 
it only for clarification when necessary. In compact ADM notation, our general action 
will look like the following: 



S - 

where S° = 1Z° 
71° 



d 4 x 



- U !JU + n$ + N£° + + S m (N, N u Nf, ...) 



(3.5) 



T° and S* = H i - V with 2 : 



7T ij 7T lm (\ 
\f9 R + — \^9ij9lr 



9il9jm 



n 



7T 



(3.6) 



The lapse and shift functions are auxiliary variables. iVj can be solved for and will 
drop out of the action while, N will enforce a constraint. With fluctuations in mind, 
we specify that: 



7T 



n ij + p ij ; N = 1 + n; $ = 0o + v 2 ; n = II + ir. 



(3.7) 



Note that we have set N = 1 in the background which we ensure by appropriate choice 
of coordinates for the reference metric, this in general implies M ^ 1. No confusion 
should arise from the common labeling ix\ = tt as opposed to the fluctuation of the 
conjugate momentum 7r of the scalar $ as what we refer to will be quite clear from the 
context. 

As mentioned, we will concentrate here on the instability of the scalar sector of 
the theory, since this is where the Higuchi bound is known to arise; to this aim, it is 



2 Note that in ADM notation y/g = and R = ^R as opposed tOy/-Wg, (4) i?. We follow [39] 
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very convenient to employ the following tensor decomposition for p iJ and for the metric 



fluctuation h^: 



1 

2 L 



3 91 J 



d l 



(3.8) 



Equipped with these conventions and notation, we now write down the most general 
expression for the massive gravity contribution, it reads: 

2m 2 [ fh 2 + S^j d 2 + (ht + h t ) 2 c 2 + e 2 (h + h t )n + f 2 N? 
-« 3 (Qii + h t ) 2 c 3 + (hf + 4 + (h + h t )ne 3 + f 3 N; 
+« 4 ((hi + h t ) 2 c 4 + (h 2 + M d 4 + (hi + h t )n e 4 + UN 2 ) ] , (3.9) 

where c n , d n , e n , / n 's are in general functions of a, b, M as of Eq. (3.1) (see Appendix A 
for an explicit expression for the functions above). 

3.1 The Higuchi bound, an example 

At this stage it seems apt to introduce the stability condition we are after, which will 
determine the Higuchi bound. In order to do so, we restrict ourselves to a very simple 
specific case encompassed by the theory above (see [38] for a detailed treatment): 

• no matter content 

• «3 = = OJ4 

• the reference metric f^ v is set equal to the background value of g^, i.e. f\ Ll/ = ^ 
z = 0, so no background contribution from the MG action 3 . 

As we will do for the general case, we write the resulting canonical action in this 
form: 

r . 1 



p' q - g p ■ K {,) ■ p ~ o q ' M(0 ' q ~ p ' V{ ' ] ' q ' 



(3.10) 



The stability conditions can be read off the matrix K and M, they should be positive 
definite for the system to be stable. As will be later explained, we will focus here 
on what is now the matrix M and disregard the so called gradient instability. In the 
simplest case, when the three conditions above are met, the Lagrangian for the helicity 
zero mode is [38]: 



Sn 



d x 



tt% + NS° + NiS* 



m 
4a 



(h ij hij - hi - 2a 2 N 2 - A^nhu) , (3.11) 



where the matter $ has been put to zero also in £°, S 1 and the MG Lagrangian starts 
at second order in perturbations. The background equations of motion simply read 
3H 2 = A, H = 0. 



3 In other words, we are projecting our theory into the region of the parameter space that gives 
back Fierz-Pauli [10]. 
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We spare the reader the simple albeit tedious algebra but list the steps required 
to reduce the Lagrangian in (3.11) to a convenient version of Eq. (3.10): 



1) Solve for jVj and plug back in the Lagrangian; 

2) Employ the field redefinitions: hij — > a 1 ! 2 hi 



p l] — > a 1 l 2 p % \ n 



-3/2 



3) Solve for the variable p l (no loss of generality); 

4) Perform the field redefinitions p l — > po + ^jjh*, h l — >• q + \h t , v 2 = m 2 — 2H 2 

5) Solve for h l by its algebraic equation; 

6) Finally, use the background e.o.m. and employ the field redefinitions below: 



p^p + H 



2H 

v 2 m 2 



(-2V 2 + 3u 2 - 3H 2 ) po 



2H / r 



u 2 m 2 



3u 2 



3H 2 )p . 
(3.12) 



The result is: 

L = Poqo - 



1 / v 2 m 2 \ 

2 \12H 2 ) 



1 {12H 2 



2 \ v 2 m 2 



Po[-V 2 



m 



9/4]p 



(3.13) 



Having switched off the scalar $ in this example, the canonical variables are just the 
couple (po,qo)', the stability condition can be read directly from the sign of the qo 
coefficient, that is: 



, 2 > <=> m 2 > 2H 2 



(3.14) 



This is the Higuchi bound [40] on m in this set up. One is of course well aware that, 
over the different cosmological epochs, it is not allowed to assume H = 0. Indeed, the 
bound as it presents itself here, serves as a warning that things need to be checked in 
more general scenarios, away from de Sitter [39]. 

We will perform steps similar to the ones mentioned above for the general theory 
but a quick route to Higuchi's result is worth it for clarifying the physical importance 
of the bound and how it comes about. Now, observationally, a lower limit on m is close 
to the last thing one worries about in MG. Quite on the contrary, the importance of 
identifying an upper bound is intuitively clear, after all, GR is known to work quite 
well. 

3.2 Higuchi VS Vainshtein 

In the m — > limit MG should reproduce GR. This does not happen quite so easily 
because the two theories have a different number of degrees of freedom. Vainshtein 
was the first to propose [13] an expansion parameter for MG solutions which makes it 
clear how inside the so called Vainshtein radius non-linearities play an important role 
in the theory. It has been proven that non-linearities help hide the additional degrees 
of freedom of MG and restore continuity with GR [55-57]. 

A heuristic way to obtain the Vainshtein radius is the following. Consider a 
generic background equation for MG, schematically, it will be of the type: 

1 



R^u + m 2 h L 



M 2 J» V 
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The linear to non-linear regime transition happens when the fluctuation ~ 1 so 
2 w„ nio„ i,™™ +w e> y2j, and ± ^ gm_ Therefore, i?„ 



that Rabcd ~ m • We also know that R a bcd 
^ m" and finally: 



rv 



M 



m 2 M 2 x 



1/3 



i'abcd 



(3.15) 



This is the celebrated Vainshtein radius. It is now clear that too large an m would 
be excluded by tests on the validity of GR in the solar system and beyond. We are 
inside ry then, and the mildest condition to impose on m would be for it to give a 
negligible contribution to the Friedmann equation of the theory. Since a cosmological 
constant contribution is allowed by data, the bound is best placed on H which for 
most of cosmic history must be consistent with GR. In other words consistency with 
cosmological observations for the period before dark energy domination requires 



2 d 
m — 

dt 



H 



H 2 



1H 3 



(3 + 3a 3 + a 4 )— - (1 + 2a 3 + a 4 )—^ + (a 3 + a^)-—^ 

tin O iln 



<£HH. (3.16) 



On the other hand we expect ^4ln(if/if ) ~ To see this, suppose that the 

reference metric describes a power-law FRW geometry. This would correspond to 
H = Ab n for some constants A and n. At the same time as we will see below we 
require 



— which implies 



b 
a 



H 

Aa r 



l/(l+n) 



(3.17) 



and so if in turn a(t) ~ t p then 



1 d 



HH/Ho 



pn + 1 
p(n + 1] 



0(1) 



(3.18) 



For a given epoch of cosmological expansion (3.16) typically amounts to a condition of 
the form 



m 



H H 2 1 H 3 

(3 + 3a 3 + a 4 )— - (1 + 2a 3 + (^4)772 + ("3 + a 4)oTH 
n n oii 



which approximately (assuming one power of H/ Hq dominates) is 

m 2 < 0(1)H 2 



(3.19) 



(3.20) 



where m 2 is the effective mass arising in the Higuchi bound (see below). It is this 
relation which creates tension with the Higuchi bound fh 2 > 2H 2 . We can of course 
imagine evading the Vainshtein bound for a specific period of cosmic evolution. For 
instance is b is chosen to behave as in a matter dominated universe, then during the 
period for which a is also matter dominated, the ratio H/Hq will remain constant, 
and so the mass term contributes nothing but an overall cosmological constant to the 
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Friedman equation, which may be absorbed into any existing contribution. However, 
without modifying b, the period of radiation domination will necessarily be modified 
from GR predictions. We will disregard such extreme tunings. 

We see now that two different checks on massive gravity, one which originates 
from studying the classical stability of the theory, the other from the necessity to make 
sense of observations and recover GR at least at early cosmic times, come together and 
force contradictory bounds on m. The order one numerical coefficients might appear 
to leave some room for maneuvering, it is therefore a good time for a more detailed 
analysis. How this tension might be resolved in ghost-free theories of MG will be the 
subject of the following sections and of a follow-up work [44, 45]. 

4 Higuchi bound: dS on dS 

In the example above, several assumptions have been made, chief among which the 
use of FP theory, one which has long been known to have instabilities, ghosts. We 
now proceed to perform an analogous analysis employing the dRGT theory of massive 
gravity in full. Let us stress again that here we are dealing with a ghost-free model at 
fully non-linear level endowed with a clear-cut and convenient hierarchy of scales. For 
the moment, we do not add matter content but make full use of the parameter space 
of the theory, including the free parameters «3, a^. What we will learn for this case 
will be qualitatively true for the most general one, i.e. the effect of the two parameters 
does not change the qualitative picture of the Higuchi- Vainshtein tension. 

We start with the Friedmann equation for the dRGT theory, without matter, 
according to Eq. (2.1): 

H 2 = \\ + m 2 {z - z 2 ) - a 3 m 2 (z 2 - z 3 /3) + \a A m 2 z 3 (4.1) 

where 1 + z — b/a. This is just a special case of the general Friedmann equation for 
an FRW reference metric 

R2 = 3Mf P + m2 ( z ~ z ^ ~ a3 m ^ z2 ~ *V 3 ) + 3 tt4 mV ' ( 4 - 2 ) 

which can easily be determined from the equation of motion (2.8). Consistency of the 
Friedmann equation and the acceleration equation a (determined from (2.8)) implies 
the relation 4 



and so 




4 We disregard as pathological the branch of solutions for which (l + 2(1 + a^)z + (a^ + a±)z 2 ) = 
[44]). 
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Putting this together the general Friedmann equation is (applicable for an arbi- 
trary FRW reference metric) 



H 1 



1 m 2 2 H 

^tjyP - (6 + 4a 3 + a 4 ) — + (3 + 3a 3 + a 4 )m — 

OlVl-n-i 6 110 



'PI 



H 2 

a 4 )m 2 j^ + (a 3 + at 4)'-^- ^g- 



m 2 if 3 
'0 



(4.5) 



We will use the Friedmann equation above in determining the Higuchi bound but 
at this stage we can already tell, because of the Vainshtein screening mechanism we 
introduced above, that the following inequality must hold before the period of dark 
energy domination: 



m 



{z - z 2 ) - a 3 {z 2 - z 3 /3) + ^a A z 3 



or more precisely 



m 



d_ 



(z - z 2 ) - a 3 (z 2 - z 3 /3) + -CK4 z 3 

3 



< HH 



(4.6) 



(4.7) 



We now proceed to give a detailed account on how to obtain the Higuchi bound in 
this case 5 . The explicit version of the calculation for the most general case will only 
be hinted at in the text of next section, while details will be reported in Appendix B. 
The explicit expression for the Lagrangian at second order in perturbations for dRGT 
without matter is: 



{dihtf m 2 hh t 3H 2 h 2 9H 2 (ht + h t ) 2 AH [ h i + f 



8a 3 



4a 



--a{pi+p t f + 2H \hipi + 



h t pt 



4a 

" « [fi + J ] 



d 2 K 



-h ijP ij - H(hi + h t )(pi + Pt ) 
2a 2 H( Pl +p t ) 



+a{hi + h t ) {H 2 -A)) 



hfA hih t A 



{ht + h t f h 2 + 



2 , hi 



4a 



2a 



A - a(hi + h t )nk 



+2 d lPi 



djhiH dih t H 



Nt + 2m z 



hi 



h 2 t 



4a 

d2 + {hi + h t ) 2 c 2 + {hi + h t )ne 2 



hi 



+f 2 N 2 - « 3 ( {hi + h t ) 2 c 3 + {hi + h t )ne 3 + ( hf + -± ) d 3 + f 3 Nf 



+a 4 (hi + h t ) 2 c 4 + {hi + h t )ne A + [ h 2 + -± d 4 + / 4 iVj 



5 There is a faster way to attain the Higuchi bound in this case, simply introducing a dressed mass 
m 2 that arises from looking at fluctuations about the de Sitter metric utilizing the de Sitter symmetry, 
and the using the representation theory statement that m 2 > 2H 2 . It is instructive and convenient 
to spell out at this stage the more general procedure, where we cannot make use of the de Sitter 
symmetry, doing it directly for the most general theory would sacrifice clarity. 



(4- 
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We refer the reader to Appendix A for the explicit expression of the functions c n , d n , e n , f n 
with the reminder that, in this case, M = 1 + z: we are in dS space. The next step is 
to solve for TVj, to obtain: 

N _ -a dj pi + dj hiH - dj h t H 
2am 2 (f 2 - h «3 + h «4) 

Plugging this back in the action, we proceed to a convenient redefinitions of variables: 

Kj -> a 1/2 hij, p ij -> a' 1/2 p ij , n -> a~ 3/2 n. (4.10) 

This will be complemented by the redefinitions on the matter scalar mode in the most 
general case. We are now ready to eliminate n, which acts as a Lagrangian multiplier, 
and use the constraint to solves for p t : 

Pt = ~Pi-§^ + ( 2 ™ 2 ( e 2 - «3e 3 + « 4 e 4 )a- 1 + (H 2 - 2A)) = 

= ~ Pl ~ mt> + {kt+ 2H )y2 ] u2 = 2m2(e2 ~ 0363 + 0464)0-1 + {R2 ~ 2A) - 

(4.11) 

It is now convenient to perform a change of variables [38]: 

v 2 h 
Pi->Po + jjjht, hi-+q + —, (4.12) 

after which one solves the algebraic equation for ht- The equation is easily solved in 
coordinate space because, as one can easily check, the coefficient of the V 2 h t term in 
the action vanishes at this stage upon using the equation of motion. Finally, in order 
to put the action in a form where there are no mixed terms of the type p ■ q (this is 
easily done in this case) we perform the canonical transformation: 

p ^ P + H(Q + A-P), qo^Q + A-P (4.13) 

where A is determined by requiring indeed that the coefficient of Q ■ P is zero. We 
refer to Appendix A for the explicit expression. At this point the action will have the 
simpler form: 

PQ-^aP' + ^Q 2 ]. (4.14) 

The last step of the analysis in then to require a > 0. Explicitly: 
m 2 (l + z)(-l + z(2 - a 3 (-2 + z) - a A z)) x 

(2H 2 + m 2 (l + z)(-l + z(2 - a 3 (-2 + z) - a A z))) > (4.15) 
There is indeed now an effective mass, 

m 2 = m 2 (l + z)(l - z(2 - a 3 (-2 + z) - a A z)). (4.16) 
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One realizes the Higuchi bound is formally identical to the one we obtained in the 
simplest case: 



~ 2 / ~ 2 

m [m 



2H 2 ) > 0. 



(4.17) 



Although this bound is easily satisfied if fh 2 < this is well-known to give instabili- 
ties in the vector sector of the theory since the kinetic term for the vector modes is 
proportional to in 2 . Consequently the real Higuchi bound is 



m 2 > 2H 2 



(4.18) 



In order to get an intuition for the Higuchi mechanism in this case, let us temporarily 
set as, «4 to zero. The form of the bound now is: 



m 2 (l - 
H 



z-2z 
H 2 

- 2 JT 2 



2H 2 > 



m 



2H 2 > 



Now it is clear what one gains when considering a reference metric which is other than 
the background value of the dynamical metric g^ v : the effective mass can be positive 
while, at the same time, the initial "bare" mass m 2 is allowed to be negative. Below 
we present the conditions that need be satisfied for the stability of the Hamiltonian. 
We split them in two convenient subsets, two branches for the value of H, assuming 
that the effective mass is positive for both . We obtain: 



r < H < 5fa 



m 



2 > 2HHj 



3H -2H 



( H > 3fi 



m 2 < 



2HH$ 



in 2 > 



(4.19) 



m 2 > 



where m 



7,2 



{3H/H - 2H 2 /H 2 ) 



m 



The branch on the left is the one which represents 



a continuous deformation of the usual case 6 , for which z = 1 <^ H = H . In such a 
branch, we recover the by now familiar m 2 > 2H 2 in the H — >■ Hq limit. This solution 
will share the same properties we elaborated upon in the first example: the bound 
itself is quite stringent and there is a strong Higuchi- Vainshtein tension (see Eq. 4.25 
below). The other branch is new: it lives in a different region of the parameter space 
where H > 3/2H , allows for a positive effective mass while keeping the bare mass 
small and negative. Let us rewrite the second condition: 



1 *\ H 




\m — 




' l H 





> 2H 2 , for H^Hn 



\m | 



> i 



(4.20) 



So, considering the freedom on the value of Ho, the Higuchi bound, by itself, is not 
that stringent anymore. It can be easily satisfied. The only caveat would be that H 



6 It is important to note here that there is no a priori reason to assume f^ u = g^ v , other than sheer 
simplicity. 
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must always be larger than Hq, which fixes the value of Hq as smaller than the current 
value of the Hubble constant. This analysis must be complemented with observations, 
i.e. MG must give very subleading contributions to the Friedmann equation. Let us 
rewrite Eq. (4.1), again setting a 3 = = a 4 , in this form: 



3|m 2 

H + 



/A + 2 |m 2 |+ Tl V 



2H (1 



(4.21) 



V J 



Clearly, we want to recover the H = a/ p/3 as m — > 0. By looking at the main 
denominator on the RHS above, one can see that the condition \m?\/HQ < 1 must be 
imposed and this conflicts with the Higuchi bound in its latest disguise as of Eq. (4.20). 

By using dRGT theory and its cherished properties, we found a new branch where 
H can live which is devoid of the typically stringent Higuchi bound. As soon as we 
focus on the Friedmann equation we realize this branch is not viable anymore. Things 
do not change qualitatively when switching on a 3 , a 4 , but we report the result below 
for completeness. Including matter is the next conceptual step, which we introduce in 
next section. The Higuchi bound with all parameters in use reads: 



m 2 (l + z)(l - z{2 - a 3 {-2 + z) - a 4 z)) > 2H 2 , (4.22) 



that is, 



2 H ( ,_ _ , , H S H 2 



m\H) = m 2 — (3 + 3a 3 + a 4 ) - 2(1 + 2a 3 + a 4 )— + (or 3 + "4)775 > 2H ■ 
-Wo V ^0 ^0/ 

(4.23) 

This of course, must be combined with the requirement due to the Vainshtein mecha- 
nism, which is obtained straightforwardly from the Friedmann equation : 

m 2 (3z - 3z 2 - 3a 3 ^ 2 + «3^ 3 + a A z 3 ) < 3H 2 (4.24) 

Combining the last 2 equations together one obtains: 

1 — z — 2z 2 — 2a 3 z — (X3Z 2 + a^z 2 + (a 3 + a 4 )r 3 



3z — 3z 2 — 3a 3 z 2 + (a 3 + a^z 



3 



> 1. (4.25) 



There exist specific values of a 3 , a^, z for which this inequality is indeed satisfied. On 
the other hand, this would require some extreme tuning on the theory while we have 
yet to consider the fact that, eventually, H and so consequently z is supposed to wildly 
vary during cosmological epochs and this alone might make seemingly special values of 
a 3 , a 4 unacceptable. It is interesting to note that both polynomials in the numerator 
and denominator are cubic unless a 3 + a 4 = in which case they are both quadratic. 
If in addition a 3 = — 1 then both polynomials are linear. Because of this fact that 
regardless of the choice of a 3 and a 4 , the polynomials entering the Higuchi bound and 
the Friedmann equation are the same order, playing with these parameters does not 
remove the Higuchi- Vainshtein tension. 
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Let us stress again that, despite using a similar analysis, the result obtained here 
differs from the ones in [38, 39] (as well as from more recent work in [41, 42]) in more 
than one, crucial, ways. First of all, the MG theory itself is a different one: the place 
of Fierz-Pauli MG is taken up by dRGT massive gravity for the reasons on which we 
expanded upon before. The second important point is that the reference metric f^ v is 
here allowed to be other than the background value of the dynamical metric g^, i.e. 
ffu, g^, the difference is parametrized by the time dependent function z. 

Before one can claim that the Higuchi-Vainshtein tension is here to stay, matter 
should be added to the picture. This is what we do below. From the scalar sector 
perspective, this corresponds to adding another scalar mode and could, in principle, 
help relax the Higuchi bound and therefore the Higuchi-Vainshtein tension. 



Consider now dRGT theory of massive gravity with FRW reference metric, and add 
matter content to it. We consider a simple expression for the matter contribution 
to the scalar sector, namely a scale field with an arbitrary potential, but we do not 
anticipate that the conceptual results will part from ours in more general scenarios. 
The idea is that our model for the scalar can mimic a generic equation of state, and 
so this should suffice for the instability analysis. 

As we have seen in the dS case, the Higuchi-Vainshtein tension forces m into small 
interval of values. It is then crucial to investigate what happens away from de Sitter, 
when matter content is included and so the theory can describe the dynamics over the 
different cosmological epochs. The full Lagrangian now reads: 



This Lagrangian has the form of (3.5). We have shown above that employing the 
freedom on the parameters as, a 4 of the theory does not change the qualitative picture 
of the stability analysis. This allows us to set 0:4 = = 03, which benefits the algebra 
as well. All of the following calculation steps may be straightforwardly generalized 
to the case where 03 and «4 are nonzero. In a future work we shall give a simpler 
derivation which is applicable in the general case and confirms the following results 



We briefly sketch the route to the Higuchi bound, which closely mimics the simpler 
one above. The main difference is that we have an additional scalar mode and therefore 
Eq. (3.10) is to be intended with matrices as oppose to bare numbers. As a consequence, 
the calculation is a bit more involved and for simplicity we do not show all the explicit 
steps in the text. We refer the reader to Appendix B for details on all the initial 
matrices entries and more. 



5 FRW on FRW 




(5.1) 



[44]. 
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The equations of motion for the full theory read: 



H 2 = m 2 (z-z 2 ) + — + - 
v ' 12 6 



7T V 



— 1 9 

7T m 



if = -— — (l — z — 2z 2 — M + 2Mz) ■ 

7T + 3Htt + V\ = ; Vi = ^L^j j vf has been defined as 7f = O ; 

z = H (M — 1 — z); a 3 = = a 4 ; • 

b H b , v 

l+Z=-=— , #0 = V77- 5.2 

a iio M o 

The algorithm one follows is similar to the one on •page 5 with a few important addi- 
tions: 

1) Solve for Ni and plug back in the Lagrangian; 

2) Employ the field redefinitions: 



hij -»■ a 1/2 hij, p lj -)■ a n a 3/2 n; 7r -»■ a 3/2 7r; <£>->■ a 3/2 y?; 

n -> a 3 ri; <9 2 /a 2 = A. (5.3) 

3) Solve for the n constraint for the variable p 1 (no loss of generality); 

4) Perform the field redefinitions: 



V 2 ut ul 1 ut V 

m h\ h ^ qo + -h, «^P1- W I 



f^Qi + ^K u 2 = m 2 (l-z-2z 2 )-2H 2 + ^-. (5.4) 

4:11 Z 

5) Solve for h l by its algebraic equation; 

6) Use the background e.o.m. and re-arrange the Lagrangian, through a series of 
canonical transformation, in such a way that it eventually takes the form: 



L = P T -Q 



aP 2 + P 2 - PQ AQo - 9l*9± + 5 Q P 1 + 9l^Ml' , 



where one can prove [39] that the second to last term can be reabsorbed 7 . The canonical 
transformations necessary to this aim are three. It is quite easy to show that they are 



7 We do not perform this step in the text because it is not necessary in order to determine the 
stability properties we arc after. 



- 15 - 



indeed canonical; it is also important to check one has enough freedom on the matrices 
that define these transformations as to be able to eventually put the Lagrangian in the 
form (5.5). 

The first transformation reads: 

p->P + fy, r (Q + a ()) -P), (5.6) 

with a, h generic symmetric 2x2 matrices. These amounts to 6 degrees of freedom 
in a, h. The second transformation is simply: Qq — > —Pq, Pq — > Qo- After these 
transformations have been applied, one rewrites the Lagrangian in the form (3.10) and 
finds the entries of the new K, M, V matrices. The requirements we impose on these 
new matrices stems from the form of Eq. (5.5): one wants the (1,1) entry of the matrix 
K to be one, the off diagonal elements to vanish and the (0,0) entry to be free of 
Laplacians 8 . That leaves three free parameters. Those can be fixed by reproducing 
the wanted value of the coefficients multiplying QiAQi, QiAQq and PiAQq, which 
saturates the number of free parameters and, at the same time, guarantees a Lagrangian 
in the form (5.5). There is a third canonical transformation, which gives its desired 
result by construction (i.e. without fixing any of the free parameters). It reads: 

P^P + A-Q, Q^Q, with A = diag[-y o/#oo,-W^ii] (5.7) 

This last step forces a final V matrix with Vio as the only non-zero entry. 

The Higuchi bound can now be read off directly from the sign of the K 00 en- 
try. Quite remarkably, upon further simplifying the expression it turns out the in- 
dependence drops out of the equations, to give a simple bound of the form: 



{3H - 2H)m 2 {3H - 2H)m 2 - 2M 2 



> (5.* 



where we have used: -g^ = H and (1 + z) = H/H . As usual requiring the vector 
modes are stable restricts this further to 

(3H - 2H)m 2 - 2HH 2 > 0. (5.9) 

This inequality immediately reduces to the usual m 2 > 2H 2 in the case without matter 
and with = g^ 9 . 

There is no need to elaborate further on the H > 3/2Hq branch at this point, as 
the discussion would be identical to the one in Section 4- As far as the H < 3/2Hq 
branch is concerned, one realizes that things are simply worsened by the realization 
that, the very same inequality m 2 > 2HHq/(3Ho — 2H) has now to hold over time. As 



8 Note that we focus on K now because the second canonical transformation flipped momenta and 
coordinates. 

9 It also gives back the result of [39] upon setting 0:3,4 to zero and using the appropriate background 
e.o.m.'s, i.e. z = 0. 
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we have seen before, making use of 03, 014 also does not change the qualitative picture. 
The general bound is (see [44] for a full derivation) 

H ( H H 2 \ 

m 2 {H) = m 2 — I (3 + 3a 3 + a 4 ) - 2(1 + 2a 3 + a 4 )— + (a 3 + 04)770 > IE 2 . 
-"0 V Mo H Q J 

and is also independent of H and the precise form for matter. It might still be in- 
teresting to describe in some detail what happens if one picks some special values for 
a 3 , 014. A nice feature is that for (a 3 + a 4 ) = the z 3 proportional contributions to 
the Friedmann equation and to the bound both vanish. At the same time, with 0:3 we 
have gained one parameter to tackle the Vainshstein screening mechanism: one can 
set the value of «3 = — 1 so that only z appears in the Friedmann equation and that 
should make things easier for large H (remember, z = H/Hq — 1). As it turns out 
after a quick calculation, the Higuchi bound, now easily satisfied upon choosing a small 
value for H , forces m 2 /H Q > H which spoils the Vainshtein back again. This pattern 
will reappear again for other tempting values of 03, 04 and one essentially notices that 
together with the freedom that comes with the a n 's , there comes a polynomial in z(t) 
which brings more conditions on the a n s in order for the Vainshstein mechanism to be 
at work. 

Although in some respects technically different from what the current literature 
offers 10 , once confronted with observations, the Higuchi bound we found still leaves 
little room if any for the value of m. 



Conclusions 

Recently, a new theory of massive gravity has been put forward [20]. Most interestingly, 
dRGT theory is free from instabilities such as the Boulware-Deser ghost and tames 
issues which have plagued previous attempts at giving mass to the graviton. The 
Higuchi bound is a condition obtained directly from the study of the classical stability 
of a given theory. It is therefore crucial to investigate how dRGT is affected by this 
bound. We have performed such a study in this paper for the dRGT model with 
an FRW reference metric, which has revealed several, interesting things. First of all, 
once stability and observations are taken into account, the Higuchi bound is still quite 
stringent on the value of m, the parameter deforming GR. In summary the bound 
states that 

H f H H 2 \ 

m 2 (H) = m 2 — (3 + 3a 3 + a 4 ) - 2(1 + 2a 3 + 04)77- + ( a 3 + "4)772 > 2H 2 . 

-"0 V -"0-^0/ 

where H is the Hubble rate of the dynamical metric and Ho that of the reference 
metric. 

This bound differs significantly from previous conditions, in particular [39] (see 
also [41-43]) in that it is independent of H and the precise form of matter. Indeed 

10 This is a different, better theory of MG, where the freedom on /„„ is also fully employed. We 
have indeed seen how there is now an additional branch which does relax the Higuchi bound. 
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the bound found in [39] m 2 > 2 (if 2 + H/N) appears to allow without constraint 
accelerating solutions. By contrast the bound found for the dRGT models is equally 
strong regardless of whether the universe is in an accelerating or decelerating regime. 
Furthermore, unlike these previous works, here we have studied a theory where the 
reference metric is distinct from the background dynamical metric. 

Generically when analyzing the Higuchi inequality, one identifies two or more 
branches of values for the Hubble parameter H that satisfy the condition. For the 
case 03 = 0^4 = 0, one branch contains the / M „ = g^ u solution in the absence of matter 
while the other one allows for a negative m 2 value whilst the effective mass m 2 , the 
one that counts, is kept positive. If one were to focus just on the instability at this 
stage, it would be correct to conclude that the branch which does not include f^ y = g^ v 
does indeed allow for a very weak, easily satisfied, Higuchi bound. However, one soon 
realizes this set of values for the parameters of the theory does not stand the test of 
observations, i.e. recovering GR wherever the Vainshtein screening mechanism is at 
work. 

Another feature that emerged in this work is that the free parameters , do 
not change the qualitative picture of the Higuchi- Vainshtein analysis; they would if one 
were frozen into a specific cosmo logical epoch, but fully considering that H = H(t) 
and that one must account for all type of matter, is enough to render the coefficients 
ineffective to this aim. 

It is interesting to note at this point how the analysis of homogeneous and isotropic 
cosmological solutions in dRGT [37] might carry the same message we learned here. In 
that work, a no-go theorem for k = FRW solutions to dRGT with Minkowski reference 
metric was presented. The way out there relied on recalling that observations only 
require approximate homogeneous and isotropic solutions, and on looking for them in 
dRGT [37, 60-64]. We plan to lift the purely homogeneous and isotropic assumptions 
on the metrics in a follow up paper [45]. As an additional note it should be said 
that the no-go theorem can be evaded by giving full dynamics to both metrics, that 
is considering multi-metric theories [46, 58, 59, 65, 66]. We did not completely rule 
out here the fact that multimetric theories could also help in relaxing the Higuchi- 
Vainshtein tension and shall return to this in subsequent work [44]. 

Verifying that dRGT theory possesses stable regions of the parameters space 
where the Higuchi- Vainshtein tension can be relaxed is quite an important task. We 
have shown that this is not the case when both the reference and dynamical metrics 
are FRW. We will report on our the results on the possible resolutions in a forthcoming 
work [44, 45]. 

Whilst this paper was being submitted, the following paper [67] appeared in the 
literature. The study in [67] very partially overlaps with the work presented here. 
The authors there focus on the de Sitter and quasi-de Sitter case (see our Section 4 ) 
for bigravity. They do not study general FRW backgrounds as we do here. Also [68] 
appeared in the literature at the same time. There the authors study the so called 
decoupling limit of solutions such as the ones described in this paper, with special 
attention to the Higuchi bound and Vainshtein mechanism in that limit. 
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6 Appendix A 

The explicit expression for the functions c n , d n , e n , f n in the MG action of Eq. (3.9) 
reads: 

3 36 3M bM\ , / 3 96 b 2 3M 36M\ 

C2= -^-7T-5--S— + ; rf 2=l — + — r - — + — - I 



("2 



3M 


bM 




+ SaJ 


8a 




6 2 \ 




2a" J ' 


h = 



4a 8a 2 8a 8a 2 J ' V 2a 8 « 2 8 « 3 4a 8q2 7 

36 2 6 3 



2(6 + aM) a(6 + aM) 



36 3M 6M\ , /l 96 6 2 3M 36M 6 2 M 



2a 8a 2 8a 4a 2 / ' \a 8a 2 4a 3 4a 4a 2 8a 

6 2 \ £ ( 36 2 26 3 
e 3 = -2a + 36 - - ; f 3 



a)' Ja V 2(6 + aM) a(6 + aM) 2a 2 (6 + aM)/' 
1 6 M bM\ ( 1 36 6 2 M 36M 6 2 M 

l 3^ _ 8^ 2_ 8^ + 8^ 2 ;' 4 ~ V 4^ + 8a^ ~ 8^ + 4^ ~ "8a^ + "8^ 

\2 + 2a)' h \2{b + aM) a(6 + aM) + 2a 2 (6 + aM) ) ' 



The explicit expression for the quantity A in Eq. (4.13) reads: 

2H (-9H 2 - 3m 2 (l + z)(-l + z(2 - a 3 (-2 + z) - a 4 z)) - 2A) 
m 2 (l + z)(-l + z{2 - a 3 (-2 + z) - a 4 z)) 

1 

X (2# 2 + m 2 (l + z)(-l + z(2 - a 3 (-2 + ij - a 4 z))) (6 ' 2) 

7 Appendix B 

We give below the entries of the matrices K, M, V as they appear before point (6) in 
Section 5, that is, before the three final canonical transformations are performed: 



o , 4ff(l+z)+2i a 4(f/+J/2+i) a 2 _J n(H+Hz+z) a 

_ | Hm 2 (l+z) 2 (-l+22) a 3f/m 2 (l+ 2 )^ 2ZA 2tf 3ff 2 (l+2)i/ 2 " 1 / 7 -p 

'S"' — \ _jf n(H+Hz+z) a m 2 7r(l-2z)(ff+ijz+z) I ' I'- 1 ; 

2tf 3H 2 (1+z)f 2 ^ 1 + 12/^2 
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5 v 2 A AHA zA 2zA 

V 00 = --H- — + 



2 2H 3H m 2 (l-z-2z 2 ) 3H 2 (1 + z) m 2 (l + z) 2 (-l + 2z) 
AHA 2 AzA 2 



V 



3m 2 u 2 (l - z - 2z 2 ) 3m 2 u 2 (l + z) 2 (-l + 2z) 
Vi 2ttA V x A ttzA v x zA 



01 



2H m 2 {l-z-2z 2 ) 3Hv 2 Hm 2 (l + z) 2 (-l + 2z) 3H 2 u 2 {l + z) 

2ttA 2 2ttzA 2 
_ 3m 2 ^ 2 (l - i - 2z 2 ) + 3Hm 2 v 2 {\ + z) 2 (-l + 2z) ( ' ' 



m 2 (l - z - 2z 2 )7i m 2 n(-l + 2z)z ttA 7fiA 
10 = ~ n 12H* + 12H~z + 3^ + 3#z/ 2 (l + i) ^ ^ 



m 2 7r\/i(l - 2z)i 3 m 2 W x (\ - z - 2z 2 ) tx 2 A tt 2 zA 

V 11 — rT o n 1" 



12# 3 z/ 2 2# 12# 2 i/ 2 6#z/ 2 6#V 2 (l + ,z) 

(7.5) 



1 



M 00 = (H {AH 2 (6tt 2 - 5m 2 (-1 + Z + 2z 2 )) + m 2 (l + z){-\ + 2z) 



x 



(-7T 2 + 2m 2 (-l + z + 2z 2 ))) + m 2 (-l + 2z) (l6# 2 - vr 2 + 2m 2 (-1 + z + 2z 2 )) z) 
2 (-6# 3 (1 + z) + #m 2 (l + 2) 2 (-l + 2z) - 3H 2 z + m 2 (-1 + z + 2z 2 ) z) A 



3Hm 2 {\ + z) 2 (-l + 2z) 

AH(H + Hz + i)A 2 



(7.6) 



3m 2 (l + 2) 2 (-l + 22)z/ 2 

m 2 Vi(-l + 2z)(H + Hz + z) n(H + Hz + z)A V X (H + Hz + z)A 



M i = Mio 



12# 3 6# 2 (1 + 3H{1 + z)z/ 



2 



7f(2F(l + z) + z)A 2n(H + Hz + z)A 2 



M 



m 2 (l + z) 2 (-l + 2z) 3m 2 (l + z) 2 (-l + 2z)v 2 

H (6H 2 {-AH 2 + vr 2 ) \/ 2 - m 2 (V? + 12H 2 V 2 ) (-1 + z + 2z 2 )) + m 2 V 2 (l - 2z)z 



ii 



12# 3 z/ 2 



A 2W l {H + Hz + z) 7r 2 (2#(l + z)+i) 

# 2 (1 + z)3 (4# 2 - 7T 2 + 2m 2 (-1 + z + 2z 2 )) 2#m 2 (l + z) 2 (-l + 2z) 
ti 2 (H + Hz + z)A 2 



3Hm 2 (l + z) 2 (-l + 2z)u 2 ' 
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Starting with generic (up to this point) entries for the symmetric matrices a and h, 
we show their final formal expression after the three canonical transformations are 
performed. We further show how to determine the various entries by the requirement 
that the final Lagrangian is of the form (5.5). First, the generic symmetric matrices: 



l ab 



hoi h u I V a 01 a u 1 



After the following canonical transformations 

Po \ ( Po \ fh 00 h i \ ( Qo\ ,( h 00 h 01 \ ( a 00 «oi \ f P o 
Vx ) ^ \Pi) \hox h n )'\Qx) \hox h n J " \a 01 J ' 

1o \ , ( Qo \ , { a oo a oi \ f Po 



Po^Qo, Qo^-Po, (7.11) 

one has to write the Lagrangian in the form (3.10) according to the new variables 
(Q, P). This, after many integrations by part, will give the new K, V, M matrices. As 
an example we report below the new K matrix. 

K o = ^oo + ^oo^oo + ft-oAi + ^oo + 2/i o(/ioi#oi + V 00 ) + 2h iV 10 , -/loo-^oi - h 01 K u - V w 
-«oi (^oo + ^oo^oo + h^Ku + M 00 + 2/i 00 (V^oi + V 00 ) + 2h 01 V w ) 

K i = -hooKoi - KiKu - V w - a i(^oo + ^oo-^oo + h^Kn + M 00 
+2hoo(hoiK i + V 00 ) + 2h 01 V 10 ) 

Kn = K u + 2aoi(h 00 Koi + h iK n + V 10 ) + ^(/ioo + ^oo^oo + h 2 01 Ku + M 00 

+2hoo(hoiK i + Voo) + 2h 01 V w ) (7.12) 

where the matrix entries K, M, V on the RHS of this last equation are to be intended 
as the "original" K,M,V entries, as given in Eq. (7.1-7.8) above. Equipped with the 
new matrices one then requires: 

- K i = 0, and solves for a i; 

- V i = and solves for hn; 

- Kn = 1 and solves for hot] 

This is clearly taking the Lagrangian towards the form in (5.5). The other conditions 
that will further fix the remaining entries (e.g. hoo = Hiiubbie, etc.) are: 

- No A's in K 00 ; 

- No A's in M i ; 

One can further massage the final expression for the Lagrangian, but, after this 
point, further requirements on a, h entries must be true by construction (i.e. the 
requirement will be a consequence of the ones above) because the freedom on the 
matrices has been saturated. We give the explicit expression for the final form of the 
K o entry in the text (Section 5) as it is by requiring it to be positive that one obtains 
the Higuchi bound. 
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